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The concept of entanglement entropy appears in multiple contexts, from 
black hole physics to quantum information theory, where it measures the en- 
tanglement of quantum states. We investigate the entanglement entropy in a 
simple model, the spin-boson model, which describes a qubit (two-level system) 
interacting with a collection of harmonic oscillators that models the environ- 
ment responsible for decoherence and dissipation. The entanglement entropy 
allows to make a precise unification between entanglement of the spin with its 
environment, decoherence, and quantum phase transitions. We derive exact 
analytical results which are confirmed by Numerical Renormalization Group 
arguments both for an ohmic and a subohmic bosonic bath. Those demonstrate 
that the entanglement entropy obeys universal scalings. We make comparisons 
with entanglement properties in the quantum Ising model and in the Dicke 
model. We also emphasize the possibility of measuring this entanglement en- 
tropy using charge qubits subject to electromagnetic noise; such measurements 
would provide an empirical proof of the existence of entanglement entropy. 
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1. INTRODUCTION 

Dissipation in quantum mechanics represents an important statistical 
mechanical problem having ramifications in systems as diverse as biological 
systems to limitations of quantum computation. A prototype model in 
this class is the Caldeira-Leggett model describing a quantum particle in a 
dissipative bath (environment) of harmonic oscillators [1]. The spin-boson 
model can be seen as a variant of the Caldeira-Leggett model where the 
quantum system is a two- level system [2,3]. Those impurity systems are 
interesting because they display both a localized (classical) and delocalized 
(quantum) phase for the spin [4,5]. This class of models were intensively 
investigated to study the emergent quantum-classical transition [6] . 
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An environment generally induces some "uncertainty" in the spin di- 
rection which is responsible for quantum decoherence, i.e., for the rapid 
vanishing of the off-diagonal elements of the spin reduced density matrix. 
On the other hand, a finite coupling between the spin and its environment 
also induces entanglement between the spin and the environmental bosons, 
i.e., the wavcfunction of the system cannot be written as a simple product 
state anymore. In this paper, we are interested in quantifying properly the 
entanglement properties between the spin and its dissipative environment. 
One good measure of entanglement for these well-defined bipartite systems 
is the von Neumann entanglement entropy [7-9]. In fact, we anticipate a 
deep connection between decoherence of the spin and (entanglement) en- 
tropy which may be seen as follows: the decoherence irreversibly converts 
the averaged or environmentally traced over spin density matrix from a 
pure state to a reduced mixture. Combining exact analytical results and 
Numerical Renormalization Group techniques applied directly on the spin- 
boson model [5,10], we are prompted to flesh out the entanglement proper- 
ties between the spin and its environment by computing the von Neumann 
entropy of the spin when integrating over the bath degrees of freedom. We 
expand our recent articles [11,12]. We will demonstrate that the resulting 
entanglement entropy constitutes an important order parameter (which ex- 
hibits universal scalings), not only of quantum phase transitions but also 
of crossovers. We will also confirm that loss of coherence of the two-level 
system (characterized either by the complete destruction of Rabi oscilla- 
tions [13] or by the strong suppression of persistent current in a ring [11,14]) 
means a prominent enhancement in the entanglement of the system; this 
enhancement of entanglement corresponds to the actual quantum phase 
transition in the case of a subohmic bath and rather corresponds to the 
dynamical incoherent crossover from damped oscillatory to overdamped 
behavior in the case of an ohmic bath [3]. For a subohmic bath, second 
order phase transitions can be interpreted as transitions where coherence is 
lost due to the emergence of a purely classical state at the transition. We 
also emphasize the possibility of measuring the entanglement entropy of the 
spin using charge qubits subject to tunable electromagnetic noise [15, 16]. 

2. SPIN-BOSON MODEL, ITS RELATIVES, ENTROPY 

In this section, we provide useful definitions. In particular, we introduce 
the spin-boson Hamiltonian as well as useful mappings towards long-range 
Ising chains [17] and Kondo models [18], and wc discuss more thoroughly 
the concept of entanglement entropy. We define properly the entanglement 
entropy E of the spin with its environment in the context of the spin-boson 
model, and wc discuss the different behaviors of E in the dclocalized and 
localized phase of the spin-boson model(s). The von Neumann entropy E 
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is strictly defined for a pure state at zero temperature; nevertheless, in Sec. 
4.3 we will briefly comment on the quantum-classical crossover at finite 
temperature where the Boltzmann-Gibbs entropy proliferates [19]. 

2.1. Spin-Boson model, Kondo physics, and Ising chains 

The Hamiltonian for the spin-boson model with a level asymmetry h 
takes the usual form [3]: 

HsB = -^(^x + ^(^z + Hose + ^0-0 X] + 

n 

where ctx and cr^ arc Pauli matrices and A is the tunneling amplitude 
between the states with cr^ = ±1. Hogc is the Hamiltonian of an infinite 
number of harmonic oscillators with frequencies which couple to the 

spin degree of freedom via the coupling constants {A„}. The heat bath is 
characterized by its spectral function: 

J{uj) = n'^XldiiJn -uj) = 2'Kaujl-'uj', (2) 

n 

where u)c is a cutoff energy^ and the dimensionless parameter a measures 
the strength of the dissipation. Similar models involving a two-level system 
coupled to a bath of harmonic oscillators are invc'stigatcd in the context 
of stimulated photon emission. A known example is the Jaynes-Cummings 
model [20]. It should be noticed that our model has a dense spectrum of 
fundamental frequencies of the harmonic oscillators. 

In the following, we will refer to the values < s < 1 as subohmic 
damping [3,6]. One standard approach is to integrate out the dissipative 
bath, leading to an effective action which is reminiscent of the classical spin 
chains with long-range correlations (here, in time): 

«5mt = j dTdT'az{T)g{T - T')a4T'), (3) 

with 

g{T) (X (4) 

at long times. For a subohmic bath, correlations in time become sufficiently 
long-range that one may expect the localization of the spin for all values of 

a > 0, at least in the limit of A ^ 0. On the other hand, by increasing A, 
one eventually expects a quantum phase transition towards a delocalized 



To simplify notations, hereafter we fix the Planck constant h = I. 
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FIG. 1. Phase diagram of the spin-boson model for the ohmic and subohmic case. 



phase for the spins. In fact, the analogy with classical spin chains leads to 
the existence of second-order quantum transitions for < s < 1 between a 
localized phase at small A and a delocalized phase at large A [17,21]; see 
the phase diagram of Fig. 1. It is relevant to mention that this transition 
has been seen by applying the Numerical Rcnormalization Group (NRG) 
method directly on the spin-boson model [6, 10, 12] but this is not visible 
in the popular non-interacting blip approximation [3]. 

In the localized regime, the tunnel splitting A between the two levels 
renormalizes to zero, i.e., the spin gets trapped in one of the states t or |, 
whereas A is strongly renormalized in the delocalized phase. 

The special case of s = 1 corresponds to the ohmic spin-boson model 
which is also equivalent to the anisotropic Kondo model; at small A, here 
the system shows a Kosterlitz-Thouless (KT) type quantum phase tran- 
sition^, separating the localized phase at a > etc from the delocalized 
phase at a < ac- To understand how the KT transition takes place as 
a function of the dissipation parameter a, one may use the analogy with 
the anisotropic Kondo model, which describes the exchange interaction 
between conduction electrons and a magnetic impurity [22] : 

kk' 

+ Y - 4iCfc'i) + hS,., (5) 

kk' 

where i?kin represents the kinetic energy of the electrons. The equivalence 
between the spin-boson model with ohmic dissipation and the Kondo model 

^In the ohmic case, the phase transition is described by RG equations similar to those 
in the XY model in two dimensions. On the other hand, the phase transition in the 
Kondo model was found earlier than the KT phase transition [22]. 
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has been first realized by Chakravarty [23] and independently by Bray and 
Moore [24]. One may identify: 

{pJi_) — and (1 + 25e/7r)2 a, (6) 

where p is the conduction electron density of states and the parameter 6^ is 
related to the phase shift caused by the Jz Kondo term and is given by [25] 

(5e = tan"i(-7rpJ;,/4). (7) 

The phase transition in the Kondo model (a purely quantum phase transi- 
tion at the absolute zero) corresponds to the (bare) Kondo coupling chang- 
ing sign. For Jz > {i.e., antiferromagnetic coupling of the Kondo spin to 
conduction electrons), the "fugacity" Jj^ is a relevant operator. Physically, 
this means that the number of spin flips proliferates, and the Kondo spin 
forms a singlet state with the conduction electrons. In the opposite < 
(ferromagnetic Kondo coupling) case, Jj^ scales to zero. No spin flips re- 
main, and the spin is "frozen" in time corresponding to the long range 
order of the Ising chain. Anderson, Yuval, and Hamman [22] succeeded in 
mapping the behavior of the Kondo model onto the solution of a long-range 
Ising spin chain. The tunneling events where the local moment spin flips 
due to scattering with conduction electrons, correspond to domain walls of 
the equivalent Ising chain. The problem is similar to ordinary tunneling of 
a two-level system. The physically new feature is the fact that conduction 
electrons represent a dissipative bath, and as result the tunneling events 
are not independent. Instead, they feature long-range interactions in time, 
and the equivalent Ising chain acquires long-range interactions. 

The equivalence between the anisotropic Kondo model and the spin- 
boson model with ohmic damping can be formulated explicitly through 
bosonization [25] . The delocalized region corresponds to the antiferromag- 
netic Kondo model, while the localized region corresponds to the ferro- 
magnetic Kondo model where the spin is frozen in time. In the delocalized 
phase of the spin-boson model, the Kondo energy scale obeys'^ [26] 

= A(A/£))"/(i-"\ (8) 

for values of a not too close to the transition and close to the transition 
Tr- assumes the exponential form of the isotropic, antiferromagnetic Kondo 
model; InTff oc 1/(q!c — o)- The critical line separating the localized and 

^In the definition of the Kondo energy scale, we have introduced a high-energy cutoff 
D for the conduction electrons which is of the order of the Fermi energy; the relation 
between the two cutoffs a)c and D will be fixed properly later; see Appendix C. 
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delocalized phase in the spin-boson model thus corresponds to ac = 1 + 
O{A/c0c) or = — I J_l|. Such results have been confirmed by applying the 
NRG approach directly to the spin-boson model [5, 10]. In the delocalized 
phase of the spin-boson model, we expect a unique nondegenerate ground 
state and the expansion in the number of spin flips in the partition function 
does not lead to a dilute gas of flips for T <C Tk which reflects the basic 
difficulty that the exact ground state is orthogonal to that for a static spin 
A = 0; consult Sec. 3.1. 



The problem of measuring entanglement in many-body systems is a lively 
field of research. Bipartite entanglement of pure states is conceptually 
well understood [9,27,28]. A useful measure of many-body entanglement 
when the total system is in a pure state is the von Neumann entanglement 
entropy [7-9]. This is obtained by focusing on bipartite systems where 
space can be divided into 2 regions, A and B. Beginning with the ground 
state pure density matrix, region B is traced over to define the reduced 
density matrix pA- From this, the von Neumann entanglement entropy 



is obtained for the subsystem A of size r. The rate at which E grows with 
the spatial extent, r, of region A is not only a fundamental measure of 
entanglement, it is also crucial to the functioning of the Density Matrix 
Renormalization Group [29]. For systems with finite correlation lengths, 
it is generally expected that E grows with the area of the boundary of 
region A [30]. In the one-dimensional case, conformally invariant systems 
have E{r) (c/3)lnr where c is the "central charge" characterizing the 
conformal field theory [31-33]; for a two-dimensional example, see Ref. 
[34]. Entanglement entropy has recently been shown to be a useful way 
of characterizing topological phases of many body theories [35 37], which 
cannot be characterized by any standard order parameter. Entanglement 
entropy is also closely related to the thermodynamic entropy of black holes 
and to the "holographic principle" relating bulk to boundary field and 
string theories [36,38]. Entanglement entropy also brings a new light on 
quantum impurity systems; a spin can be entangled to a bath of conduction 
electrons [39] or to a dissipative bosonic environment [11,12,40,41]. For a 
quantum impurity system (A is the spin): 



2.2. Entanglement entropy 



E{r) = -Tr[pAlog2 pa], 



(9) 



E = -P+ l0g2 P+ - P- Iog2 P- , 



(10) 



where 




(11) 
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For the spin-boson model, {oy) = because Hsb is invariant under Oy 
—ay. Moreover, for the spin-boson model, we get the equalities: 



where £g is the ground state energy of Hsb- Since Hsb and Harm 
are related by a unitary transformation, they have the same ground state 
energy (up to an unimportant constant) . The field h couples directly to the 
spin in both models, so we have {S^) = {a ^12. A similar relationship does 
not hold between {S^ and (da:), and thus E does not (exactly) measure 
the entanglement between the Kondo impurity and the conduction band. 

At a = 0, the qubit is decoupled from the environment; thus, {gx) = 
A/V/i^ + A^ and (a^) = -/i/V/i^ + A^. With p+ = 1 and p- = 0, we 
check that = for all values of A and h. Moreover, at large h, we also 
must have E ^ because the qubit is localized in the state with (cr^) = — 1 
and (ax) =0. In a similar way, in the localized phase, since the qubit is also 
localized in one classical state at arbitrarily small h then E 0^. Below, 
we want to understand how E interpolates between all these limits. We 
argue that the entanglement entropy of the spin with its environment will 
allow us to make a unification between entanglement of the spin with its 
environment, decoherence, and quantum phase transitions. Entanglement 
near a quantum phase transition is of great current interest [9,42,43]. 

It is also relevant to note the relation between entanglement and quantum 
decoherence of the spin. Quantum decoherence implies a rapid reduction 
of the off-diagonal elements of the spin density matrix, i.e., {<Tx) = which 
results in p± 1/2 and i? ^ 1 if (cr^) vanishes at small h (delocalizcd 
phase). We will show below that for the ohmic case, the quantum decoher- 
ence or maximal entanglement occurs at the Toulouse limit a = 1/2 [44] 
whereas for the subohmic case this rather happens at the phase transition. 



After these generalities, we consider the ohmic situation and a large 
level asymmetry A//i <C 1 where perturbation theory applies and a corre- 
spondence to the well-known P{E) theory for dissipative systems can be 
formulated [45]. We are also interested in extracting the spin observables. 




-2d£g/dA 
2d£g/dh, 



(12) 
(13) 



3. OHMIC CASE: SPIN OBSERVABLES AND 
ENTANGLEMENT 



3.1. Large level asymmetry and P(E) theory 



*This result is intuitive if one absorbs the coupling with the environment into the 
tunneling term [14]. In the localized phase, the tunneling term is renormalized (almost) 
to zero which then ensures a disentanglement (decoupling) with the environment. 
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FIG. 2. Two-level system with a large level asymmetry A/h <g; 1. 

For A = 0, the spin is completely localized in the | state. Now, for finite 
A, we want to evaluate the probability pdei that this spin down electron 
flips up. In the absence of the environment, it is straightforward to obtain 
that pdei = A^//i^ + ... . Below, we will incude the effect of the environment 
and see how this quantity gets modified. It should be noted that pdei plays 
an important role because it is related to {<Tz), through 

(a,) + 1 - 0{pdei). (14) 

For A = 0, one gets two classical states and the Hamiltonian for either 
of the states depends on the bosons: 

n n 

We can easily handle this Hamiltonian when A„ = 0. The stationary 
wavefunctions are those with a fixed number of bosons n„ = ajja„ in each 
mode corresponding to the energy X^n'^n"-"- 7^ 0, one may absorb 

the linear term in the redefinition 

an,^ = a„ + ^ (16) 

In 

an,l — fln ^• 

where 7„ = A„/a;„. Since the bosons have been shifted, the vacua for the 
two states are not the same. More precisely, let us consider the vacuum 

a™,i|Oi)=0. (17) 

This is equivalent to 



(18) 
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We then observe that |0|) is not the vaccum for the bosons hnked to the 
state t- This is rather a coherent state and thus 

\0,) = exp (-^) g (2^|{n„},). (19) 

For convenience of notations, in this expression we have assiimed that m is 
fixed. To first order in A, this impHes that the ground state |0|) acquires 
corrections proportional to all possible states |{»^m}T)- 

\9) = \Oi) + ^M\{nmh), (20) 

with 

^{n„} = aIM!^. (21) 

h + rimWm 

The probability that the spin flips up, then takes the form: 

Pdel = \^{nm]? + - ■ (22) 

By carefully summing over all modes m, one can make a formal link with 
the P{E) theory [45] of dissipative tunneling problems (see Appendix A): 

If there is no dissipation, then P{E) = S(E), and thus this reduces to the 
known qubit result p^e; — A^/ft,^ + ... . In the case of Nyquist noise or ohmic 
dissipation where P{E) oc E'^"~^ [45,46], the delocalization probability 
rather evolves as 

A V /i ^ 



Note that Pdei is small at high h ^ lOc and one reproduces Pdei ~ A^/h^. 
On the other hand, for a < 1, one observes that pdei increases for smaller h 
and eventually reaches the maximum value pdei 1- This shows that the 
problem becomes highly non-perturbative at small h ior a < 1; one indeed 
expects a delocalized phase for a < 1 and a localized phase for a > 1. It 
is important to note that for small level asymmetries, the point a = is 
unstable, i.e., perturbation theory in A already breaks down at small a. 
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The result in Eq. (24) is in agreement with the expansion of the ground 
state energy to second order in A (we consider the case where h> 0) [14] : 



where T is the incomplete gamma function. In the regime of relatively large 
h (but still small compared to the high-energy cutoff uic), one may simplify 



r 1 - 2a, — ~ r(l - 2a) - 



{h/ujc) 



l-2a 



1 - 2a 



This immediately leads to: 



2 ■ 4 Vwc 
1 



(26) 



(27) 



1 - 2a 



r(l - 2a) 



2q-1 



Thus, {uz) = 2d£g/dh can easily be derived. At relatively large values of 
the level asymmetry, one gets: 



1 /A 



2 Vw, 



(1 - 2a)r(l - 2a) 



2q-2 



(28) 



We check that (cr^) + 1 ^ 0{pdei)- In a similar way, one can extract 
(fx) = —2£g/dA, resulting in: 



A 



1 



- r(l - 2a) 



(29) 



l-2a 

It should be noted that the limit a —> 1/2 is always well-defined, since 

(30) 



iim(l-r(.),-^ 



1- 1-2/"" 1 
hm = iny. 



In particular, this leads to: 

{o'x)a=l/2 = 



X— >0 X 



■A in 



(31) 
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in agreement with the exact results obtained from the non-interacting 
resonant level model at the Toulouse point a = 1/2 [44]; see Appendix B. 
This implies the important equality between the two cutoffs of the theory: 

D{a = 1/2) = — . (32) 



4. RESULTS FOR A VANISHING LEVEL ASYMMETRY 

In the dclocalizcd phase, one can use the Bcthc Ansatz solution of the 
anisotropic Kondo model or of the equivalent interacting resonant level 
model; consult Ref. [47,48] and Appendix C. While general expressions are 
quite complicated [14], it is instructive to derive simple scaling iorinis for the 
observables in the limits h <C Tk and h ^ Tk [11]. We will also apply the 
NRG method which shows an excellent agreement with the Bethe Ansatz 
results; NRG parameters have been defined in our Ref. [12]. 

4.1. Spin observables 

For h 3> Tk, wc check that the Bethc Ansatz calculations reproduce the 
perturbative results of Sec. 3.1. For the sake of clarity, technical details 
are given (hidden) in Appendix C. For h <C Tk, we obtain, 

/ N 2e^(T^ r[l + 1/(2- 2a)] / /i \ 

.^^^-^^ = - r[l + a/(2-2.)] [t-k) ' 

where 

6 = alna-|-(l-a)ln(l-a). (34) 

Note that {uz) oc H/Tk at small h, in keeping with the Kondo Fermi liquid 

ground state. The local susceptibility of the spin converges to 1/A for 
a ^ 0^ in accordance with the two-level description and diverges in the 
vicinity of the KT phase transition as a result of the exponential suppression 
of the Kondo energy scale. It is relevant to observe that the longitudinal 
spin magnetization {gz) only depends on the "fixed point" properties, i.e, 
this is a universal function of H/Tk in the delocalized phase. Exactly at 
the KT transition, one can adapt the results of Anderson and Yuval [49] 
and the result is that (az) jumps by the amount y^l/oc along the quantum 
critical line where ac = l + 0{A/u)c)- This results in a non-universal jump 
~ 1 — 0{A/ojc) at the KT transition [5]. Remember that in contrast to the 
universal jump of the superfluid density in two-dimensional XY models, in 



SPor a ^ 0, r[l] = 1, r[3/2] = V7r/2, exp(6/(2(l - a)) = 1, and Tk = A. 
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FIG. 3. Longitudinal spin magnetization (cz) as a function of the detuning h from 
Bethe Ansatz techniques (top) and NRG calculations (bottom). 



the spin-boson model, the jump in the longitudinal magnetization is non- 
universal for finite A/cJc- Finally, far in the loealized phase, one rather 
predicts (az) « — 1 + ©((A/wc)^) [40]. In Fig. 3 we present two curves of 
(az) as a function of h, which have been obtained with the Bethe Ansatz 
and with the NRG calculations, respectively. One can notice the very 
good agreement between the two approaches. By increasing the detuning 
h, the abrupt jump in (az) occurring at the quantum phase transition is 
progressively replaced by a smooth behavior; at finite h there is a smooth 
crossover separating the delocalized and the localized regime. 

The leading behavior of (ax) in the delocalized phase takes the form: 

hm (a.) = + C,{a)^, (35) 



with 
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FIG. 4. The transverse spin magnetization {ctx) versus a for different level asym- 
metries h obtained from the Bethe Ansatz calculations (top) and from NRG (bottom). 



,-V(2-2o) r[l-l/(2-2a)] 
^^("^^ V^(l-a)r[l-«/(2-2a)]- ^^^^ 

As a — > 0, Tk —f A and Ci(0) — 1, so we recover the exact result 
(o'a;)a=/i=o = 1^- As we turn on the coupUng to the environment, we 
introduce some uncertainty in the spin direction and {a^) decreases. It 
should he noted that {a^) does not only depend on fixed point properties; in 
the delocalized phase, (ax) still contains a perturbative part in A/lOc! 

For a < 1/2, the monotonic decrease of Tk/ A dominates. In contrast, 
for a > 1/2, the first term in Eq. (35) dominates and we have 

{<7x)a>l/2,h^0 ^ 7^ 7 • (37) 

2a — 1 Wc 



"Exactly at o = 0, the (small) first term in Eq. (35) is not present; see Appendix C. 
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In fact, this result can also be recovered using the perturbation theory 
of Sec. 3.1 (consult Eq. (29) for — > 0). This result clearly emphasizes 
that the observable (ct,.) is continuous and smMlaX the KT transition in the 
ohmic spin-boson model. This is also consistent with the work by Anderson 
and Yuval which predicts {gx) A/wc exactly at the phase transition [49]. 

The fact that {a^) becomes very small ^ A/wc in the localized phase can 
be understood from the mapping between the anisotropic Kondo model and 
the Ising model with long-range correlations in time. Here, the partition 
function of the system can be expanded in powers of [3,22,49]: 



A2" r^drr f^dT2n T-r 



,2n 



where the function W describes the (long-range) interaction between kinks 
(spin flips) located at positions and tj, and A = A/uic- 

The scaling procedure requires to modify the short-time cutoff Tc = 
u)~^ Tc — dvc- Two important features then appear by modifying Tc- 

the dependence of W ^ can be included in a global renormalization 

of A and configurations with a kink-antikink pair at distance between Tc 
and Tc — dTc have to be replaced by configurations where this pair is absent. 
The number of removed pairs is proportional to dTc and the (free) energy 
variation is roughly — A'^ {tc)t~'^ dTc'. 

^^=-A\tc)t-\ (39) 

This is equivalent to change u)c ^ i^c + dA = A: 

BP 

^=A^(A). (40) 
The renormalization of A in the ohmic case is [3] : 

This leads to 

F(A) = / dAA^A). (42) 

The ground state energy can he recovered by taking the limit A ^ 0. This 
is continuous at the quantum phase transition. Extending this result to 
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the delocalized phase, we find that for a > 1/2 this expression is efi'ectively 
convergent at low energies, and essentially one recovers the result obtained 
by Bethe Ansatz or by perturbation theory: 



This ensures that in the dclocahzcd phase, for a > 1/2, (a^) is of the order 
of A/cOc- This result is definitely emphasized in Figs. 4. In fact, it should 
be noted that (a^) is always continuous at a quantum phase transition. 

Finally, we can also check that (cTj,) evolves continuously close to a = 
1/2. In the limit a 1/2, one can take Ci(a) = (4/7r)r(l - 2a) 
4/(7r(l — 2a)) and use the identity D{a = 1/2) = Auc/tt to find {ax) — > 
— (4/7r)-\/rff/-Dln(Tx/-D), in agreement with the "non-interacting" reso- 
nant level description valid at the very specific point a = 1/2. 



Now, we are able to study thoroughly the entanglement between the spin 
and its environment, at least in the context of an ohmic bath; see Fig. 5. 

First, in the delocalized phase, a.t h = 0, E is entirely determined by {(ix)- 
Thus, at small a, E grows linearly with a and for a > 1/2, (a^) « A/dj^ ~ 
which ensures that E(a > 1/2. h = 0) — > 1, i.e., p± ~ 1/2; see Fig. 5. It 
should be remembered that in the ohmic case, the reduction of quantum 
coherence (superposition) of the spin rather occurs at a = 1/2 and not 
at the quantum phase transition. Additionally, the maximal entanglement 
occurring at a = 1/2 allows to explain the beginning of the "incoherent" 
dynamical behavior [3]. More precisely, as a < 1/2, one has damped 
Rabi oscillations at low temperatures; the Rabi oscillation frequency A and 
damping rate F are given by A = cos riTx {a) and F = sin tjTk {a) with [3] 
T] = 7ra/(2(l — a)). In contrast, for 1/2 < a < 1 the Rabi oscillations 
disappear and the behavior becomes completely incoherent [13]. When 
the entanglement with the environment becomes too important one cannot 
distinguish in which state, t or 4, is the two-level system anymore! 

The plateau of entanglement is reminiscent of the maximal entanglement 
between the spin and its cloud of electrons in the SU(2) Kondo model [9,39]. 

Second, in the localized phase, we obtain (ctz) « —1 and (ax) « — 
and therefore i? « — for infinitesimal h = 0^ . Since dissipation localizes 
the spin in the "down" state for /i ~ 0+, we do not expect the entropy to 
depend strongly on the external field; the perturbation theory in A predicts 
E ^ — (A/a)c)^ ln(A/cL'c) to leading order for all h [40]. As wc approach 
the phase transition from the localized side, this behavior is replaced by 
E ~ -(A/cJc) In(A/a,V;) in agreement with (ct^) = -l + C(A/wc) at the KT 
transition. This implies a non-univc;rsal jump of the entanglement entropy 
at the KT phase transition. This non-universal jump is reminiscent of the 




(43) 



4.2. Entanglement entropy 
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non-universal jump in the longitudinal magnetization (a^). It should be 
noted that at finite h, the sharp non-analyticity at the quantum phase 
transition is replaced by a maximum which signifies the crossover regime 
h ~ Tk [11]- Finally, we like to emphasize that in the delocalized phase, the 
Kondo scale controls the entanglement between the spin and the bath. For 
h ^ Tk, the disentanglement proceeds as [h/TKY whereas for h ^ Tr-, 
E vanishes as (Tftr/Zi)^"^", up to logarithmic corrections. This universal 
behavior stems from (tr^) since {ax) evolves very smoothly with the level 
asymmetry h] see Appendix C. Although E is defined at zero temperature, 
this universality is reminiscent of thermodynamic quantities [50,51]. 
To be more precise, for a small level asymmetry, we find the scaling, 

lim E[a,^,h)^E{a,\Q)-h{a)(J^\ , (44) 
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FIG. 6. Coefficients fci and k2 as a function of a. 



where the coefRcient fci (a) is given by, 

2eT^ / r[l + l/(2-2a)] y 
TT\n2 \r[l + a/i2^2a)]) ' 



and 



S(a, A,0) = 1 - 



For a large level asymmetry this rather results in 



lim E{a,A,h) = fcafa) ( — 



In 



Tk 



where the pre-factor is given by 

(l-Q)e-''r(3/2-Q) 



V^ln2 r(l-a) 



0.2 0.4 0.6 0.8 1 

a 



(45) 



(46) 



(47) 



(48) 



The fact that ^2 ^ when a ~ 1 is consistent with the fact that E does not 
depend much on the level asymmetry in the localized phase. The evolutions 
of the coefficients ki and ^2 as a function of a are shown in Fig. 6. At the 
Toulouse point [44], one can easily check that 



fci(a = l/2) = ^— - 
TT^ m 2 



(49) 



fc2(a = 1/2) 



7rln2' 
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in agreement with the non-interacting resonant level model of Appendix 
B. In Fig. 5, the full solutions of {ax) and (ct^) obtained with Bethe Ansatz 
and with NRG have been used to plot E versus a for different h. 

4.3. Delocalized phase and Boltzmann-Gibbs entropy 

There is a prominent entanglement entropy in the delocalized phase. 
Even though this von Neumann entanglement entropy is generally not a 
good measure of entanglement at finite temperature [9] , it is certainly rele- 
vant to discuss the evolution of the entropy of the spin at finite temperature. 

For temperatures much smaller than Tk, the system is almost in a "pure" 
state since all the excited states only have a small weight ~ cxp —{Tk/T}. 
Thus, one does not expect a dramatic modification in the entanglement 
entropy of the system for temperatures smaller than Tk ■ This can be seen 
by considering the entropy of formation of the system [9] . Entanglement at 
finite temperature has been evidenced in different materials [52]. On the 
other hand, if the temperature becomes much larger than the Kondo energy 
scale, the entropy of the spin should be of purely "thermal" origin [53] and 
should converge to^ S = — Tr[plog2p] (normalized to 1) where 



here, Ha is the Hamiltonian for the spin alone and Za the associated par- 
tition function. The main coupling with the reservoir becomes of thermal 
origin. Now. it should be noted that just above Tk the coupling between 
the bath and the spin produces some corrections in the thermodynamic 
properties of the interacting resonant level model which produces a reduc- 
tion of the thermal entropy of the spin (above Tk) [50,51]: 



In this sense, the Kondo temperature corresponds to the energy scale at 
which the thermal entropy of the spin is strongly reduced and is replaced by 
a prominent entanglement entropy between the spin and its environment. 
It should be noted that the thermal entropy S{T) exhibits a similar power- 
law behavior as the von Neumann entropy E{h) at zero temperature. 

In the localized phase, there is no entanglement entropy, and thus the 
spin is only governed by its thermal Boltzmann-Gibbs entropy. 



'^Tlic P(E) theory allows to prove this rigorously. For large temperatures, the (small) 
coupling a does not aflect much the spin obscrvables since the P(E) function converges 
to that of the uncoupled two-level system (Appendix D). Zero-point fluctuations of the 
bath become clearly negligible in front of thermal excitations of the bath oscillators. 




(50) 




(51) 
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5. SUBOHMIC CASE: CRITICALLY ENTANGLED SYSTEM 

Below, we rather focus on the subohmic situation < s < 1 which 
exhibits a second-order quantum phase transition [6, 10, 12, 54] by analogy 
to classical spin chains [17,21] (as a function of A). The case s — 1/2 
is of particular interest since it can be realized through a charge qubit 
(dot) subject to the electromagnetic noise of an RLC transmission line [55] ; 
the ohmic case corresponds to an LC transmission line [4, 14]. In the 
limit of small A, the P{E) theory of Appendix A suggests that the spin is 
localized for all a. On the other hand, a quantum critical point Ac (a) [6,12] 
still exists and a delocalizcd phase for the spin occurs for A > Ac where 
Ac ^ when a ^ 0. It should be mentioned that, at small s, some critical 
exponents are different from those in the classical Ising model [6,56]. 

The order parameter for this dclocalizcd-localizcd or quantum-classical 
transition is usually the longitudinal magnetization of the spin (cr^) [6] 
which now vanishes continuously at the phase transition, 

((7.) (X (Ac - A)^, (52) 

where simple hyperscaling relations easily lead to (see Appendix D): 



2/3 = - s), (53) 

where the correlation length exponent u obeys [21] — a/2(1 — s) for 
s ^ 1 whereas at small s, one finds [6] l/i' — s through a small s-expansion. 
For the subohmic spin-boson model, f >2, which ensures that /3 > 0. 

To better characterize those second-order quantum phase transitions, 
we are prompted to examine the entropy of entanglement shared between 
the spin and its environment. The subohmic model cannot be solved ex- 
actly. Then, we will resort to the "bosonic" NRG [5,6,10,12,57] and to 
hyperscaling relations which are applicable if the fixed point is not trivially 
Gaussian [58]. Technical details are given in Appendix D. In fact, there is 
a simple way to conceive that the maximum of E should coincide with the 
quantum phase transition for the subohmic situation. Starting from the 
delocalized phase, since the longitudinal magnetization (ct^) = at ft = 0, 



OA 21n2 ^ 



1 



< 0, (54) 



where wc have introduced the transverse susceptibility x± = d{(Tx) /dA; 
this expression is still valid at Ac since for the subohmic case (ctz) is zero 
at the quantum critical point. Since x± = d{(Jx)/dA and {<Tx) are positive 
quantities, this implies that dE/dA < in the delocalizcd phase. In the 
localized phase, E is rather controlled by the finite longitudinal magneti- 
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FIG. 7. Susceptibility Xz(^) in the subohmic spin-boson model. 



zation and by the susceptibility Xz - 
dE _ X.IK) I 



5A 21n2(CT^ 



In 



9|((T^)|/aA > (see Fig. 7)*^: 
'\ + {o 



1 - (a,) 



> 0. 



(55) 



Xz is related to the rapid localization of the spin. Thus, dE/dA > in the 
localized phase. Eqs. (54) and (55) imply that the entanglement entropy 
is maximum at the phase transition. This shows that those second-order 
impurity quantum phase transitions are always accompanied by a maximum 
(cusp) in the entropy of entanglement. In general, such zero-temperature 
impurity critical points show a "fractional" entanglement entropy which 
depends on the dissipation strength a through (ax) and thus is not uni- 
versal. In fact, starting from the localized phase, by analogy to the ohmic 
case, we find (ax) ~ c(Q;)A/a-'c where c((x) increases by decreasing a. 

On the other hand, E unambiguously exhibits universal scalings even 
though the entanglement is two-sided, so that two numbers are necessary to 
specify E (one for each of the two ways of approaching Ac). More precisely, 
near the qantum critical point Ac, the transverse spin susceptibility obeys 



X±(A) -X±(Ac 



lA 



(56) 



where the exponent C is defined as C = v — 2. For the subohmic spin- 
boson model, one finds C > for all < s < 1, ensuring that x_l does 
not diverge at the transition. Taking into account that i^Ux) is continuous 
at the transition, Eq. (56) thus implies that E always rises linearly for 



'Here, we neglect the less relevant contribution from xi 
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FIG. 8. NRG results for the case s = 1/2, which might be realized through an RLC 
transmission line. In the subohmic case, E shows a cusp at the phase transition. The 
NRG reproduces the universal exponent k = 2/3 ~ 0.66 for different values of Ac. 



A — > — means that we approach the quantum critical point from 
the delocahzed (locahzed) region. It should be noted that the coefhcients 
C-i- and c_ can be different in the delocalized and in the localized phase 
and especially at s = 1/2 [12]. Through the NRG, we have also checked 
that c+ < 0, emphasizing that in the delocalized phase XJ- substantially 
increases at A+. For all < s < 1, this strongly underlines the duality 
between the enhancement of entanglement and the strong reduction of the 
two-spin state quantum superposition near the phase transition. 
In the localized phase, we obtain the scaling behavior: 



X.(A)cx |A-A,ri+'^(i-^)/2+a; 



(57) 
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here a ^ when v{l — s)/2 > 1, and we identify a = Xz{^c)- ^'^^ ^ > 
Xz = 0. Using a small s-expansion [6] which predicts v = 1/s B.i small s, 
we deduce that Xz diverges at A~ for s > 1/3, which is well verified with 
the NRG approach; sec Fig. 7. Using Eqs. (52) and (53), then^ 

i?(A,)-£;(A)oc |A-A,r(i-^). (58) 

We observe that the decay of the von Neumann entropy E in the localized 
phase is faster than linear for all s > 1/2 and the behavior becomes strictly 
linear at s = 1/2. The entanglement entropy at s = 1/2 is shown in Fig. 
8. In the limit s ^ 1, we also check that E becomes rapidly siipprcssed at 
A~ which is a reminiscence of the KT transition (ohmic case, s = 1). 

Now, we shall discuss the scaling of E with the longitudinal field. Using 
the hyperscaling relations, for all < s < 1 we find that (see Appendix D): 

£(/i,Ae)-i;(0,Ae)oc-|/ir, (59) 

and 

« = I = Ktt^)- <»' 

where S is defined in a usual way as: {(Tz){^, Ac) cx It is interesting 

to observe that one can always expand p± (e) = p±{h = Q)± mh?^^ at small 
h and m > to satisfy dhE{h) < (the field h favors a product state). It 
is also relevant to notice that in the delocalized phase, E decreases as h"^ 
similar to the ohmic case, whereas in the localized phase by approaching 
the phase transition we rather find a linear decrease of E with h. 

For the subohmic case, the entanglement entropy allows us to establish 
important connections between impurity entanglement, quantum decoher- 
ence (or strong reduction of the quantum superposition of the two spin 
states) when approaching the phase transition from the delocalized phase, 
and rapid disentanglement in the localized or classical phase for the spin 
(the spin is rapidly frozen in one classical state due to dissipation) . 

We note that the entanglement entropy E of the spin is an interesting 
new order parameter of those second-order (impurity) quantum phase tran- 
sitions even though the maximum of entanglement at the quantum critical 
point is not universal (but settlings ttre universal). 

6. COMPARISON WITH OTHER MODELS 

Here, we compare the entanglement entropy of the subohmic spin-boson 
model, which yields a parallel with long-range Ising spin chains and ex- 

^For s < 1/3, one rather finds E{Ac) - E{A) oc | A - Ac| + 5 . 
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hibits a second-order phase transition, with the entanglement entropy of 
the quantum Ising model [42] and with the one of the Dicke model [9,59]. 

6.1. One-dimensional Quantum Ising model 

In fact, we argue that the cusp occurring in the subohmic situation of 
the spin-boson model is very typical of second-order phase transitions and 
of quantum Ising models. Remember that the spin-boson model can be 
mapped onto an effective action where the spin is subject to long-range 
correlations I/t^"*"** in time. Interchanging space and time, one thus ex- 
pects that the the single-spin (single-site) entanglement properties in the 
quantum Ising model yields a very similar cusp at the phase transition. 

As an example, let us consider the one-dimensional (ID) quantum Ising 
model which is described by the following Hamiltonian: 

H=-E{r^^Ui+^j)' (61) 

where (7° is the ath Pauli matrix (a=x,y or z) at site j, and A describes the 
Ising coupling between spins whereas the transverse field A has been set to 
1. The model can be solved exactly using the Jordan- Wigner transforma- 
tion [54] and the system exhibits a second-order quantum phase transition 
at Ac = 1 separating a paramagnetic phase where (c7|) = from a ferro- 
magnetic phase. In the ferromagnetic phase, the order parameter obeys 
(we use the fact that the dynamic critical exponent z = 1 and v = 1 [54]): 

(<7|)cx(A-l)V8. (62) 

Here, /3 = 1/8 as a reminiscence of the classical two-dimensional Ising 
model [54]. By analogy to the spin-boson model, one can verify that (crj) 
is continuous at the second-order quantum phase transition [42]: 

/XV 1 r , 1 + ACOS0 

(o-?) = - / d(P , ^ . (63) 

^ TT Jo ^1 + y + 2Xcos(l> 

In the delocalized phase (A <C 1), we recover that (aj) 1 whereas in 
the localized phase (A » 1), we observe that (frj) is small ~ 1/A. Similar 
to the spin-boson model, one can check that the transverse susceptibility 
X_L = d{aj)/dX does not diverge at the transition, which implies that in 
the delocalized phase, the single-spin (single-site) entanglement entropy 
rises linearly close to the quantum critical point [42]. In the localized 
region, from Eq. (55), we also check that dE/dX oc (A — 1)5"^ diverges at 
Ac, meaning that E goes very rapidly to zero in the localized phase. To 
siimmarizc!, the ID quantum Ising model yields a single-site entanglement 
entropy which is very similar to that of the subohmic spin-boson model. 
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6.2. Dicke model 

Now, we briefly discuss the entanglement properties [59] in the one-mode 
superradiancc (Dickc) model [60] where collective and coherent behavior of 
pseudospins (atoms) is induced by coupling — with interaction A — to a 
physically distinct single-boson subsystem. The Hamiltonian reads [61,62]: 

H = uioJz+(^a^a+^=(a + a^){J+ + J-), (64) 

where this form follows from the introduction of collective spin operators 
of length j ~ N/2, in the thermodynamic limit {N,j) oo, the system 
undergoes a rjuantum phase transition at a critical coupling Ac = y/Luujo/2, 
at which point the system changes from a large unexcited normal phase to 
a superradiant one in which both the field and atomic collection acquire 
macroscopic occupations [61]. In the thermodynamic limit, the problem 
reduces to a two-mode problem by using the Holstein-Primakoff transfor- 
mation [63] of the angular momentum operators Jz = {b^b — j), J+ = 
b'^sj^j — b'^b, and J_ = here, b and b^ are bosonic operators. In the 
normal phase (A < Ac), expanding the square roots directly this explicitly 
results in: 

= LJob^b - jujo +uja^a + X{a + a^){b + b'<). (65) 

The problem is solvable going to the position-momentum representation 
for the two oscillators, x = (l/\/2a7)(a + a^) and y = (l/^/2a^)(6 + b^), 
with the momenta defined canonically. After diagonalizing the problem, 
one gets two independent (eff'ective) oscillators, and the energies e^?^ of the 
two independent oscillator modes read [61]: 

(4"^) ' = ^ + ± V('^2_^2)2 + 16AW) . (66) 

Crucially, one can see that e^-* remains real only for A < Ac; thus ff^") 
remains valid only in the normal phase. In this phase, the ground state 
energy is given by —jujo which is ^ 0{j) whereas the excitation energies 
e^'* are 0(1), implying that the excitation spectrum is quasicontinuous. 
It should be noted that this normal phase has a definite parity, since the 
Hamiltonian commutes with the parity operator 11 = cxp(?7r[a^n+6^fe]) [61]. 

In the superradiant phase (A > Ac), the field and the atomic ensemble 
acquires macroscopic occupations, thus one has to redefine: 



a'^ + V^, b'^ - v^. 



(67) 
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where a and /3 are 0{j). One can still diagonalize the problem; the solution 
of the two-independent oscillator modes reveal the excitation energies [61]: 




+ 4a;2a;2, (68) 



where 



2 



M=^ = -. (69) 

Again, one can check that the excitation e'l'' remains real as long as A > Ac, 
and the ground state energy is given by —j[{2X'^/u)) + {u)gLo/8X'^)]. It should 
be noted that the global symmetry 11 is broken at the phase transition. 

In fact, one may identify e_ with the excitation energy of a photon branch 
and e_|_ with the excitation energy of an atomic banch: e_ vanishes at the 
superradiant transition whereas e+ remains finite. Moreover, for A — > Ac, 
from either direction one may identify: 



where v = 1/4 is the critical exponent describing the divergence of the 

— 1/2 

characteristic length ^ = e_ ' . The fact that e_ vanishes at Ac implies 
that the Dicke model also exhibits a second-order quantum phase transition. 

In the case of an infinite number of atoms N oo, the entanglement 
entropy between the atoms and field follows the critical behavior [59] (the 
reduced density matrix is computed from the ground state wavefunction): 

£;<x, = -i^log2|A-Ac| =log2^, z/=l/4. (71) 

The entanglement between the atoms and field diverges with the same 
critical exponent as the characteristic length — a clear demonstration of 
critical entanglement. For finite N , Et^t at the critical point varies as log2 N 
by analogy to conformal field theories in (1 + 1) dimensions [31]. 

The normal phase allows a ferromagnetic ordering for the pseudospins 
{{Jz) — > — j), whereas in the superradiant phase, (J^) decreases progres- 
sively and continuously. It might be worth to revisit entanglement proper- 
ties for an infinite number of oscillators, e.g., by investigating entanglement 
properties of the dissipative quantum Ising chain [64-66] ; in this case, the 
spin order parameter vanishes at the quantum phase transition. 
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7. CONCLUSION AND MEASUREMENT 

A dissipative environment induces a mixed state for the spin (some un- 
certainty in the spin direction) at zero temperature, which is responsible 
for quantum dccohcrcncc and entanglement entropy. Wc have shed light 
on the deep connection between entanglement entropy of the spin with its 
environment, quantum decoherence, and quantum phase transitions. 

7.1. Brief summary 

An ohmic bosonic environment allows to perform exact calculations by 
resorting to the Bcthc Ansatz solution of the interacting resonant level 
system [47, 48] . Recently, an effort has been done to generalize the Bcthe 
Ansatz approach for interacting resonant levels out of equilibrium [67] . For 
a weak coupling with the environment, the eigenvalues of the spin reduced 
density matrix obey p+ = 1 — 0{a) and p_ = 0{q) [27] which ensures 
that the entanglement entropy of the spin with its environment rises lin- 
early with the coupling a. On the other hand, at a = 1/2, the off-diagonal 
elements of the spin reduced density matrix become very small ^ A/wg 
and as a result the system exhibits a plateau at maximal entanglement 
in the delocalized phase. For the ohmic case, the quantum decoherence 
of the spin does not occur at the KT phase transition but rather at the 
Toulouse limit [44] which also marks the beginning of the incoherent dy- 
namical behavior where Rabi oscillations disappear [3, 13]. In contrast to 
the longitudinal spin magnetization which is only controlled by the nature 
of the Kondo (Fermi-liquid) fixed point, the transverse spin magnetization 
also depends on irrelevant operators which characterize the localized phase 
(creation of a single kink-antikink pair in time; see Sec. 4.1). At the KT 
phase transition, the entanglement entropy exhibits a non-universal jump 
which is reminiscent of the non-universal jump in the longitudinal spin 
magnetization (cr^) [5,40]. In contrast, {ax) is continuous at the phase 
transition. In the localized state, the spin lies in one classical state and 
almost disentangles from its environment; p+ — > 1 and p_ — > 0. In the 
delocalized phase, we have also proven that the Kondo energy scale con- 
trols the entanglement properties of the spin-boson model; for h <C Tk the 
entropy decreases as {h/TKY accordance with the Fermi liquid ground 
state (or with (ctz) cx H/Tk) whereas for h » Tk, the entanglement goes 
slowly to zero as {Tk /hY~'^'^ . In fact, the entanglement entropy shows 
universal scalings that is reminiscent of thermodynamic quantities [50,51]. 

In the subohmic situation, the system exhibits a second-order quantum 
phase transition (as a function of A) and we have shown that the entan- 
glement between the spin and its environment is always enhanced near 
the quantum critical point. More precisely, the entanglement entropy of 
the spin always exhibits a visible "cusp" at the quantum phase transi- 
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tion and the quantum decoherence of the spin is progressively destroyed 
in the vicinity of the phase transition. In this case, we have established 
clear connections between impurity entanglement, strong reduction of the 
quantum superposition of the two spin states when approaching the phase 
transition from the delocalized phase, and rapid disentanglement in the lo- 
calized (classical) phase for the spin. The entanglement entropy represents 
an important new order parameter of those second-order quantum phase 
transitions even though the maximum of entanglement entropy at the quan- 
tum critical point is not universal (but, scalings arc universal). A similar 
cusp occurs in the single-site entropy in the quantum Ising model [9,42], 
emphasing the link between spin-boson and Ising models [17,21]. 

For a finite level asymmetry h, the entanglement entropy always yields 
a maximum characterizing the departure from the delocalized phase. 

Open questions concern the time-evolution of entanglement and the link 
between entanglement entropy and entropy production out of equilibrium 
in quantum dot systems [68]. It would also be worth to study the entan- 
glement of a chain of quantum spins with a bosonic bath characterized by 
a dense spectrum [64-66] , and to compare with the Dicke model where the 
entanglement between the (pseudo)spins and field diverges at the (second- 
order) phase transition as E = log2 ^ where ^ is the characteristic length. 

7.2. Measurement 

The problem of measuring entanglement in many-body quantum systems 
represents an active field of research [27,28]. An important open question 
in the study of quantum entanglement is whether it can be measured ex- 
perimentally. In fact, the model considered here is realized in noisy charge 
qubits, composed either of a large metallic dot [4] (the single electron box) 
or a superconducting island [15] (the Cooper pair box). The gate voltage 
controls the level asymmetry h, and A corresponds to the tunneling am- 
plitude between the dot and the lead or the Josephson coupling energy of 
the junction. If the gate voltage source is placed in series with an external 
impedance, voltage fluctuations will give rise to dissipation even at zero 
temperature [4] . An external resistor, which can be modelled by a long LC 
transmission line, would describe the ohmic situation (s = 1) whereas an 
RLC transmission line could mimic the subohmic case where s = 1/2 [55]. 
The parameter a can be varied in situ when a two-dimensional electron gas 
acts as the ohmic dissipative environment [16]. In the spin-boson model, 
E depends only on {(Tx) and {(Jz), so it can be constructed from physi- 
cal observables. While these quantities would obviously be measured at 
finite temperature, it is possible to recover the ground state density matrix 
by extrapolating them to their zero-temperature values. Charge measure- 
ments [69] yield the quantity {(Tz), which represents the occupation of the 
dot or island. In a ring geometry, the application of a magnetic flux gener- 
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ates a persistent current that is proportional to {ux) [14]. Another promis- 
ing system is the atomic quantum dot, which also permits experimental 
control of the coupling between the dot and the bosonic reservoir [70]. 
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APPENDIX A 

Spin-boson model and P(E) theory 

Here, we formulate a precise correspondence between the perturbative 
calculation at small A/ft, <C 1 and the well-known P{E) theory of dissipa- 
tive tunneling problems. From Sec. 3.1, by definition P{E) obeys: 

P(iJ) = ne-l^"'l Y,^-^^^^5{E-nmnu:^). (A.l) 

m rim 

We have explicitly taken into account that P{E) is contributed by all the 
modes m. This expression can be rewritten as 

Pi^E) = /'^]J^b:^^e-l^"l'e^(^-""'^'")* (A.2) 

m rim 



17m I _^-inmOJmt-\'ym\'^ 



m rim 



|2/-„-iu;m*_ 



J 27r 



^ iEt Kit) 
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where we have defined K{t) = J{t) — J(0) and 



/J, , \ 2 

^^^<5(a;-a;™)e--* (A.3) 



Jo ^ 
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At long imaginary times r, in the ohmic case, this gives rise to the usual 
function [45,46]: 

/'+00 7 

K{t) = -2a _(i-e-'^l^l)e-'^/'^= (A.4) 
Jo ^ 
= -2aln{l +Uc\t\). 

The logarithm appears due to the summation over the infinite number 
of modes. At low energy E, this results in P{E) oc E^"~^; there is an 
orthogonality catastrophe for a < 1 whereas for a > 1, P{E) — > at 
^ 0. P{E) can be interpreted as the probability to "emit" the energy E 
to the bath when flipping the spin. At finite temperature, one rather gets: 

K{t) = a dco^{n{\uj\) + ^(a;))(e-''^* - 1). (A.5) 

J — oo ^ 

The step-function 9{lv) corresponds to the zero-point fluctuations and n{u>) 
is the Bose-Einstein distribution function. If the temperature is large, K{t) 
is controlled by thermal excitations of the bath oscillators. The real part 
of K{t) thus becomes ~ —aTt in the long-time limit. The bath function 
P{E) becomes a Lorentzian and for very small a the result converges to 
that of an uncoupled two-level system. 

In the subohmic case, P{E) is more rapidly suppressed at low E which 
suggests that the spin is always in the localized phase at small A. 



APPENDIX B 

Non-interacting resonant level at a = 1/2 

For a = 1/2, the spin-boson model is known to be equivalent to a non- 
interacting resonant level model. The level Green function obeys: 

Gdi^n) = -. j—if^ r^, (B.i) 

iu)n-h + iTKsgn{u)n) 
and we have introduced the Kondo energy scale at a — 1/2, 

= -p-. (B.2) 
The density of states of the impurity thus takes the form: 
f \ ^c^ r< f \ ^ TKsgn{uj) 

TT TT (W — hy + J-x 
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The associated ground-state energy thus reads: 



dio 



(B.4) 



The cutoff D is important to regularize the integral. Now, we can perform 
the change of variable u ^ u + h which results in: 



D 



2tt w2 



. T2 



duj 



'-K 



D 



. T2 



h 
2' 



(B.5) 



Finally, 



h 

— tan 



-h 



(B.6) 



Thus, this results in: 



{cTz) = — tan ^ 

TT 



-h 



sgn{h) 



(B.7) 



= tan 

TT 



Tk' 



In a similar way, we identify 



2 /TV, /Tf+^ 




(Bi 



It should be noted that at a= 1/2 results from perturbation theory agree 
with the non-interacting resonant level picture if we equate D{a = 1/2) = 
AlOc/tt; this equality will be proven in Appendix C. 
First, in the limit h <^ Tk, we find 



[CTz 



-{2/n){h/TK) 



(B.9) 



-(4/7r)v/TV7D[l + HTk/D)]. 



The result for (ctz) is consistent with the Kondo ground state, where S is 
fully screened and {Sz) oc H/Tk at small h. Since both (ax) and (az) are 
small, the system is close to maximal entanglement, and assuming Tk "C D, 
we find the following scaling behavior of the entanglement entropy: 



lim £;(l/2, A,/i) = £;(l/2, A,0) 



h<^TK 



7r2ln2 \Tk 



7^] , (B.IO) 
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where 



1 + ln 



TV 
D 



(B.ll) 



The second term in Eq. (B.IO) is a universal function of h/Tx, with a 
quadratic dependence on energy that arises from the Kondo Fermi hquid 
behavior of {<Tz)- In the opposite Umit h T^, we find that 

{a,) ^ -l + 2TK/{nh) (B.12) 
^ -{4/n)^/T^ln{h/D). 

Again, it is instructive to observe that the leading /i-dependence of E has 
a universal form dictated by (ctx): 

lim E(l/2, A, h) = — ^ In ( —] . (B.13) 

The two scaling functions obtained for h ^ and for h ^ imply 
that the Kondo energy scale controls the entanglement between the 
spin and its environment. On the other hand, because of the logarithmic 
correction, E approaches zero slowly at large h. 



APPENDIX C 
Bethe Ansatz approach 
Here, we use the Bethe- Ansatz solution of Ponomarenko [48] on the in- 
teracting resonant level model which was inspired by the rcsonant-lcvcl 
model constructed by Filyov and Wiegmann [47] (even though a differ- 
ence emerges in the number of channels and in the method of (5-function 
regularization). In fact, (a-^) has been calculated by Ponomarenko: 

_^Y-oo (-1)" ^^^^ r[l + (2n+l)/(2--2a)] / h A ^"+^ , ^ 
I \ Z^n=0 (2n-H)n! V[\+a{2-a+\) I {2-1a)\ \Tk j ' ^ 1 \ 

— ^ , .^„,„,„ ^ , / . \-1n{\-ci) \^-^) 



with the parameter b = a In a -I- (1 — Q!)ln(l — a). In the limiting case 
h -C Tk, this results in (the term with n = dominates): 

lini la ) - 2e^ r[l + 1/(2 - 2a)] 
hK^J"'' - T[l + a/{2 - 2a)] [tk ) ' > 

Those expressions are consistent with those at a = 1/2: (cr^) — > — (2/7r)(/i/TR: ) 
for h-^Tx and {a^} -1 + 2TK/{Trh) at /i > Tk- 
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{ox) can be obtained by differentiating the ground state energy £g with 
respect to A. Since {cfz) = 2d£g/dh, the ground state energy is given by 

_ Eo+Ei{h),h<TK 

Eo + Ei{h = TK)+E2{h),h>TK, 

where Ei (h) and E2 {h) are obtained by integrating the appropriate pieces 
of (cr^): 

Ei{h) = i t dh'{<j,)h'<T., E2{h) = i /" dh'{a,)h'>T^. (C.4) 
^ Jo ^ JTk 

We can determine the integration constant £^0 by comparing the ground 
state energy to the perturbative result a,t Tk- 

First, we need to evaluate the functions Ei{h) and E2{h). We obtain: 

1 „ (—1)" (2n + l)b 

Em = (2. + 2)(2i+l)n! -" (^•^) 

r[l + (2n+l)/(2-2a)] ^"+2 
r[l + a(2n+l)/(2-2a)] V^^if, 

In particular, 

^,(, = T.) = -^ (2n + 2)(2i+l)./ "^ (^-^^ 

r[l + (2/; + l)/(2-2a)] ^ 
r[l + a(2n+l)/(2-2Q!)] ^' 

The second integral gives 

pm- ^-^^ 1 V°° (-1)" -„b r[l/2 + (l-a)n] 



1 - 2n(l - a) 



^ X l-2n(l-a) 



In the limit h T^, equating the coefficients of the h?°' ^ terms obtained 
with perturbation theory and Bethe Ansatz calculations, one identifies 

2a)^ - ^ , r(3/2~a) 11 

This matches Eq. (26) of the paper by Ccxiraschi and Biittikcr [14]. More- 
over, the energy Eq is obtained by comparing the ground state energy at 
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large h to the perturbative result. For a ^ we precisely get^: 
F-^'^/'^V 1 1 r[l-l/(2-2a)] ^ 

The first term is a remnant of the localized phase. 
Gathering all these informations, for h < Tk, we find: 

(a.) = ;r^- + Ci(a)^ + IkfJ±^a), (C.IO) 



where 

Ci(a) 



2a- 1 ' A A ' \Tk 



e 2(T^ r[l - 1/(2 - 2a)] 



(2„+l)b 

1 ^ (-l)"eW^ 



This general (and quite complicated) form of (a^) is consistent with the 
computation of the persistent current by Cedraschi and Biittiker [14]. 

Now, we will seek to simplify this expression assuming that h ^ Tk- 
{ax) evolves very smoothly ~ {h/Tx)^ with the longitudinal field (level 
asymmetry). Additionally, it should be noted that in the limit /i — > 0^: 

(-.)(c.^O) = ^^+C;.(a)^. (C.ll) 

This shows that in the delocalized phase, (ctx) is not only a universal func- 
tion depending on the fixed point properties, i.e., on Tk, but it still contains 
the perturbative term in A/cOc typical of the localized phase. In fact, the 
second term becomes negligible for a > 1/2 which ensures that for a > 1/2: 

(^x) = ^-^^. (C.12) 
2a — 1 Wc 



-"^It should be noted that for a = 0, we can't set e^l'^'^ = 1 in the expression (25); 
instead we must take e^/'^" = 1 + K/llJc- This produces a term oc {h/ujc)^" which 
becomes a constant at a = and exactly cancels the term in Eq. 

^For a ^ 0, the function Ci 1, which also ensures that {(Jx) — ► 1 (and as discussed 
above, the first term in Eq. (C.ll) or in Eq should not be present when a = 0). 
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Moreover, the limit a ^ 1/2 must be taken carefully: Ci(l/2) = (4/7r)r(l— 
2a) — > 4/(7r(l — 2a)), so the two terms combine to give 

{a^)^-{4.|^:)^/T^\n{TK/D), (C.13) 

where we have used D{a = 1/2) = Aujc/t^- One can check that this agrees 
with the solution at a = 1/2 to leading order. 
In the opposite limit, h > T^, we find 



where 



- a) n!(l - 2n(l - a)) 



X 



r(i + (l-a)n) 
r(l " an) 

y2 



^ ^^ri\{l-2n{l- a)) T{l-an) 

1 



X ( 2nyi-2"(i-«) ^ 

1 — a 



For h ^ Tif , we can again simplify. 



to_^M = -l+(l^)c.(a)(2l)"° (0^15) 



where 



2e-^ r(3/2-a) 
V^(l-2a) r(l-a) • ^^-^^^ 



Using the precious relation between the two high-energy cutoffs D and ujc, 
one can check the good agreement with the perturbative results of Sec. 3.1. 
When a 1/2, {(ix) contains two additional terms which conspire with 
the other terms to produce the logarithm obtained by perturbation theory; 
we do not write them because they cancel each other for a ^ 1/2. 
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APPENDIX D 
Subohmic case and Second-order phase transitions 

For a second-order impurity quantum phase transition we can apply the 
foUowing scahng ansatz for the impurity part of the free energy, 

Fi„p=T/(|A-Ae|/TV-,/iT-''), (D.l) 

to relate the critical exponents associated with E to other critical exponents 
such as the correlation length exponent ;y. Even though the temperature T 
is introduced for the scaling analysis, the entanglement entropy E (which 
is defined for a pure state) is evaluated at zero temperature. The crossover 
from the quantum critical regime to one or other of the stable regimes, 
defines the energy scale h* that vanishes at as: 

h* oc |Ae - A]*"^, (D.2) 

It should be noted that the ansatz (D.l) is usually well justified when the 

fixed point is interacting [6, 12,58]; for a Gaussian fixed point the scaling 
function would also depend upon dangerously irrelevant variables. 

Integrating out the boson degrees of freedom induces a long-range in- 
teraction in time which results in the following term in the action [6], 
<5mt = / drdr' az{T)g{T — T')az{T'), where (?(r) oc I/t^"*"** at long times. 
Since the action is dimensionless, this implies that cr^ behaves as ~ r'*"^^/^. 
Moreover, from the term which contains the level asymmetry we also infer 
that h behaves as ~ t'^-'^-''-)/'^ or r(^+'')/^. Thus, this results in: 

6=^. (D.3) 

Defining the exponent S as {az){h,Ac) oc one also finds: 

<- 1 + s 



1 



(D.4) 



This is consistent with the NRG results which predict that the local sus- 
ceptibihty Xz — d\{az)\/dh at the quantum critical point diverges as T~*. 
In the small s limit, this leads to 1/5 k, 1 — 2s + 0(,s^); this can also be 
recovered by resorting to a small s expansion [6]. This shows that the 
transverse field term is not important to characterize critical exponents. 

Another important exponent is (3 which (here) is defined as (ctz) oc (Ac — 
A)'^ in the localized phase. Using the analysis above this results in: 

= yhS-^ = y f • (D.5) 
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(c7^) = -2dF/dA obeys A^) - (a^)(0, A^) (x -\h\^/^ with 

^/S=Y^JI-1M. (D.6) 

In fact, since l/S{s) > 2/6{s), we can check that the scahng of the en- 
tanglement entropy with h mainly stems from the behavior of the longitu- 
dinal spin magnetization. On the other hand, it is interesting to observe 
that, since at small s the critical exponent v obeys v = 1/s, one also gets 
1/6 = 2/5 = K, which is verified through the NRG for < s < 1/2 [12]. 

All critical exponents obey hypcrscaling which has been verified by NRG 
[6, 10, 12, 57]. We can thus safely conclude that the critical fixed point is 
interacting for all < s < 1. Moreover, the exponents /3 and 6 are distinct 
from the ones of the long-range classical one-dimensional Ising model (in 
coordinate space) with l/r^+* interaction for s < 1/2 which displays mean- 
field behavior there. The naive quantum-classical mapping fails for the 
subohmic spin-boson model if interactions are long-range in time [6,56]. 

Now, let us define the susceptibility Xz = —d\{<Jz)\/dA in the localized 
phase; we then obtain: 

X.(A)oc |A-A,|-i+''(i-«)/2 + a; (D.7) 

here a ^ when i/(l — s)/2 > 1, and we identify a = Xz{A~). Since v = 1/s 
at small s, we infer that Xz(A) should diverge at Ac for s > 1/3. This is 
well verified through NRG; see Fig. 7. Moreover, near the quantum critical 
point Ac, the transverse spin susceptibility x± = 29^F/(9A)^ obeys 

X±(A) = x±(Ac)+c+/_|A-Ac|f, (D.8) 

and from the free energy defined in Eq. (D.l): 

C = u-2. (D.9) 

For the subohmic spin-boson model, one finds i/ > 2 for all < s < 1, 
ensuring that x± does not diverge at the transition. Taking into account 
that (a^) is continuous at the transition, Eq. (D.8) thus implies that E 
always rises linearly for A A+ (A^ means that we approach the quan- 
tum critical point from the delocalized (localized) region) [12]. It should be 
noted that the coefficients c_|_ and c_ can be different in the delocalized and 
in the localized phase. Through the NRG, we also check that c+ < 0, em- 
phasizing that in the delocalized phase x± substantially increases at A+, 
and that x± shows a clear singularity for 2 < u < 3 or 1/3 < s < 0.94 [12]. 
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